A random fiber bundle with many discontinuities in the threshold distribution 
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We study the breakdown of a random fiber bundle model (RFBM) with n-discontinuities in the 
threshold distribution using the global load sharing scheme. In other words, n + 1 different classes 
of fibers identified on the basis of their threshold strengths are mixed such that the strengths of 
the fibers in the i — th class are uniformly distributed between the values (T2i-2 and a2i~i where 
1 < i < n + 1. Moreover, there is a gap in the threshold distribution between i — th and i -\-l — th 
class. We show that although the critical stress depends on the parameter values of the system, 
the critical exponents are identical to that obtained in the recursive dynamics of a RFBM with a 
uniform distribution and global load sharing. The avalanche size distribution (ASD), on the other 
hand, shows a non-universal, non-power law behavior for smaller values of avalanche sizes which 
becomes prominent only when a critical distribution is approached. We establish that the behavior 
of the avalanche size distribution for an arbitrary n is qualitatively similar to a RFBM with a single 
discontinuity in the threshold distribution (n = 1), especially when the density and the range of 
threshold values of fibers belonging to strongest (n -\- l)-th class is kept identical in all the cases. 



PACS numbers: 46.50. +a, 62.20.Mk, 64.60.Ht, 81.05.Ni 



I. INTRODUCTION 

What causes fracture of materials in nature? Are 
there any precursor that signals the imminence of a com- 
plete breakdown so that we can avoid them taking place? 
These are few of the many questions that physicists and 
engineers are looking at now a days to explore the frac- 
ture dynamics of heterogeneous materials [l|, l^]. A prior 
knowledge of failure properties of such materials are of 
extreme importance for problems related to physics of 
breakdown, material science as well as architectural, me- 
chanical and textile engineering. The simplest of all the 
attempts is a model of fibers with randomly distributed 
threshold strengths known as the random fiber bundle 
model. 

In a random fiber bundle model (RFBM) d, 0, 1, 
i, 0, S i, 0, [H E, E Q m, H [i3> fibers with 
stochastically distributed values of threshold strengths 
(i.e., the maximum external stress a fiber can withstand) 
are clamped at both the ends. Since the threshold of a 
fiber depends crucially upon the presence of defects in 
that particular fiber, it is indeed useful to assign ran- 
dom threshold strength to each fiber of the bundle. The 
threshold strengths however are chosen from a given dis- 
tribution, usually approximated by a Weibull or uniform 
distribution. Under the application of a weak external 
load, the fibers with threshold lying below the applied 
stress, break and the resulting additional load is dis- 
tributed among the intact fibers using a load sharing 
rule. This redistribution of stress causes further failures 
and the dynamics stops when the system reaches a fixed 
point at which no further failure takes place. To resume 



the recursive dynamics, the external load is further in- 
creased to break the next weakest intact fiber. The pro- 
cess continues until the complete breakdown of the entire 
bundle. 

In this work, we shall use the global load sharing 
(GLS) rule where the additional load generated due to 
the breaking of a fiber is shared equally by all the intact 
fibers of the bundle. In the GLS scheme, the breakdown 
of the fiber bundle can be interpreted as a continuous 
phase transition with well defined critical stress and crit- 
ical exponents 0, 0, i, [l3|- For a RFBM with GLS, it 
has also been established that the distribution -D(A) of 
an avalanche of size A, defined as the number of fibers 
broken between two successive loadings, satisfies a uni- 
versal power-law ^(A) ~ A~^/^ in the limit of A ^ cx3 
J^ 3| ■ In a recent work, Pradhan, Hansen and Hemmer 
ll( showed that in the vicinity of the critical distribution 
where the average external load on the bundle is maxi- 
mum, the avalanche size distribution shows a crossover 
to a new power-law behavior given by £'(A) ~ A"'^/^ for 
small values of A << Ac while for A >> Ac, the A^^/^ 
behavior is recovered. The characteristic size Ac around 
which the crossover occurs diverges in a power-law fash- 
ion as the critical distribution is approached and at the 
critical distribution the A"'^/^ behavior is observed for 
the entire range of A. 

In a recent paper, the authors 15] investigated the ro- 
bustness of the above universal power-law behavior of the 
avalanche size distribution by introducing a discontinuity 
in the threshold distribution p{(Jth) which is given by 
1 



p{(Jth) = 



((72 -O-l) 



1 



< ath < CTi 

O"! < (Tth < <72 
0-2 < O-th < 1- 



(1) 
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1 - ((72 - CTi 

Hence, two types of fibers (weaker and stronger) sepa- 




FIG. 1: Mixed Uniform Distribution with {n + 1) classes of 
fibers and n-gaps or discontinuities. 



rated by a gap in their threshold distributions coexist in 
the same bundle with more than half of the fibers belong- 
ing to the stronger class. It has been established that 
there exists a non-universal non-power law behavior in 
the avalanche size distribution for small A which crosses 
over to the universal behavior A~^'^ in the asymptotic 
limit of A. Most interestingly, it was pointed out that the 
non-universality becomes prominent only when a critical 
distribution is approached. The threshold distribution 
given in Eq. 1 with CT2 = 0.5 is the critical distribution 
of a RFBM with single discontinuity [l5| . because as soon 
as the redistributed stress reaches 0.5, the bundle is crit- 
ical and breaks down completely with an infinitesimal 
increase in the external load. At a critical distribution, 
however, there is a crossover from the non-universal be- 
havior to a power-law behavior with A~'^/^ for large A. A 
recent study of an infinite gap generalization of the above 
model predicts a new exponent (=9/4) of the avalanche 
size distribution [17|. 

The natural question that remains is what would 
be the effect of many such discontinuities on the 
avalanche size distribution? We address this issue in the 
present communication where we study a RFBM with 
ri-discontinuities in the threshold distribution and inves- 
tigate its effect on the critical behavior and the avalanche 
size distribution. The paper is organized as follows: In 
Section II, we introduce the model and derive the crit- 
ical stress and exponents using the recursive dynamics 
approach . The results on the avalanche size distribution 
and comparison with the n = 1 case studied previously 
are presented in section III. We make concluding com- 
ments in section IV. 



II. THE MODEL AND THE RECURSIVE 
DYNAMICS 

The threshold distribution of a RFBM with n- 
discontinuities studied in this paper is shown in Fig. 1. 
The threshold range of the weakest and the strongest 
class of fibers is from to cti and (T2n to 1, respectively. 



The mathematical form of the normalized threshold dis- 
tribution is given by 



Picrth) = 



1 



1 -i]r=i('^2i -(J2i-i) 

= :; ^F^ for (72,-2 < CTth < CT2j-i(2) 

= otherwise 

(3) 

where 1 < i < n+1, gi = cr2i-o'2i-i and cto = 0, a2n+i = 
1. The above distribution suggests that (n + l) different 
classes of fibers with the threshold of the fibers of the i-th 
class ranging from cr2i-2 to (T2i-i are mixed. At the same 
time the restriction < cti < (T2... < a2i-i < <J2i < 1 
ensures the existence of a gap or discontinuity given by 
gi in the threshold strengths between i-th and {i + l)-th 
class of fibers. We also assume that the threshold val- 
ues of the i-th class of fibers are uniformly distributed 
within the range <72i-2 to <72i-i for all z, a condition that 
leads to a set of (n -I- 1) additional relations connecting 
the parameters of the system in the following way. Dis- 
tributing a fraction ft of total number of fibers to the 
i-th class (with ^" fi — ^) and using the uniformity 
condition mentioned above we get 



/.- 



f^2i 



cr2i 



1 - Er=i 9t 



(4) 



such that (To — 0, a2n+i = 1 and once again 1 < i < ^+ 1 
The model therefore involves {2n) values of CTi's, {n + 1)- 
values of density fi, the conditions given by Eq. 4 along 
with the additional restrictions ^" fi — 1. Therefore, 
the total number of free variables that can be chosen 
independently reduces to 2n. It should also be noted 
that since 1 — X]i=i 5* i^ always less than unity, we must 
have ((721-1 - 0-2,-2) < fi for aU i. 

We shall estimate the critical stress of the above model 
within the framework of a recursive dynamics and global 
load sharing. The fibers belonging to all {n + l)-classes 
cooperatively participate in sharing the additional load 
arising due to the breaking of the weaker fibers. The 
fraction of unbroken fibers after a time step t+1, (ienoted 
by Ut+i, is related to Ut through the relation |10l.[l5l| 



Ut+i = 1 - P{<Jt) = 1 






(5) 



where P{<Jt) is the fraction of broken fibers with the ap- 
plied stress (7 and redistributed stress at, and is given 
as 

rcrt 

P{crt) = / p{(Jth)dath- 
Jo 

Similarly, the redistributed stress after a time step {t+ 1) 
satisfy the recursive relation 



Ct+l 



Ut 



t+i 



a-Pi'^t))- 



(6) 



The fixed point solution for U {= U*) at which no further 
failure takes place can be obtained by solving the above 
recursive relations (5) and (6)[i3]- Assuming that the 
redistributed stress at some instant t, at exceeds a2n (i-e, 
when the redistributed stress initiates the breaking of the 
{n + l)-th class fibers), we obtain 



U^ 



t+i 



l-[ 



I 

4 = 1 



1 



1 



SO that we get the fixed point solution 
1 ,, 



-(- 



U* = 



1- 



(T2n)] 



(7) 



2(1 -Er 5.)^ 

Along the same line of arguments, we find the redis- 
tributed stress at the fixed point 



* 1 1 

a = 

2 2 

Using Esq. 
given by 



1 . 



(8) 



(7) and (8) we find that the critical stress is 



1 



(9) 



^ 4[i-(Er50]- 

which is the applied load per fiber at which half of the 
fibers break. When the external load exceeds CTc even 
by an infinitesimal amount, there is no real solution of 
U* and a* which signals the complete break down of the 
bundle. 

The critical stress of the mixed model varies with the 
gaps Qi and in the limit g^ — > for each i, we retrieve 
the critical stress Cc = 1/4 for a RFBM with uniform 
distribution [iO|. We also get back the result of single 
discontinuity case i.e., n = 1 [1^ if only gi ^ 0. The 
equation (8) also shows that the redistributed stress at- 
tains the maximum value of 0.5 at the critical external 
stress (Tc. 

Calling the redistributed stress as x from now on- 
wards, we consider the constitutive equation, F(x) — 
Nx{l — P{x) where F{x) is the average external load 
when the redistributed stress is x. The load F{x) max- 
imizes when the redistributed stress x is equal to 0.5. 
Since the maximum value of the redistributed stress is 
equal to 0.5, (T2„ must be less than 0.5 so that some fibers 
from {n + l)-th class also fail at the critical point. The 
constraint equations (4) suggest that for recursive dy- 
namics to hold good, more than half of the fibers must be- 
long to the (n-l- l)-th class with thresholds lying between 
a2n to 1. We therefore define the threshold distribu- 
tion p{(Jth) (Eq. 2), along with the constraint conditions 
and a2n = 0.5 as the critical distribution, ft can also be 
shown in a straightforward way that the order parameter 
exponent (/3) and the susceptibility exponents (7) stick 
to the mean field (GLS) values with /3 = 7 = i/2, even 
in the presence of an arbitrary number of discontinuities 
provided uniformity condition is satisfied and hence the 
critical behavior remains unaltered although there is an 
appreciable change in the critical stress. 



III. AVALANCHE SIZE DISTRIBUTION (ASD) 

We shall now turn our attention towards the ASD 
of RFBM in the presence of many discontinuities in 
the threshold distribution which is the key point of our 
study. We show below that the discontinuities have a 
non-trivial affect on the ASD when a critical distribu- 
tion is approached. We also argue that a situation with 
many discontinuities is qualitatively similar to the single- 
discontinuity case. The scenario is established consider- 
ing a special case with n = 2, where gi = (J2 ^ o'l and 
<?2 = 174 — CT3 so that three classes of fibers with range of 
thresholds lying between (0 to txi), ((T2 to as) and {a^ to 
1), respectively, coexist in the bundle while /i, /2 and /a 
are the corresponding densities. 

Below is shown some of the allowed distributions which 
satisfy the restrictions mentioned in Eq. 4 and case 4 
refers to a critical distribution. 



Case 


/i 


/2 


CTi 


(^2 


c^s 


0-4 


a-c 


1 


0.10 


.20 


0.08 


0.16 


0.32 


0.44 


0.31 


2 


0.15 


0.25 


0.135 


0.15 


0.375 


0.46 


0.27 


3 


0.05 


0.05 


0.04 


0.16 


0.20 


0.28 


0.31 


4 


0.10 


0.20 


0.07 


0.16 


0.30 


0.50 


0.35 


5 


0.10 


0.2 


0.08 


0.26 


0.42 


0.44 


0.31 


6 


0.20 


0.i5 


0.16 


0.2 


0.32 


0.48 


0.31 



We shall now generalize the results of Hemmer and 
Hansen [J| to study the avalanche behavior of the n = 
2-discontinuity model. The general expression for the 
avalanche size distribution with GLS is given as 



N 



.A-l 



A! 



dxp{x 



(1 - «(^)) 
a(x) 



X exp({-a(a;) -H lna(a;)}A) (10) 



where x is the redistributed stress and the upper limit 
of the integration {xc) is the redistributed stress at the 
critical point. Also, a{x) = xp[x)/l — P{x) is the number 
of fibers that break as a result of breaking a fiber with 
threshold strength x by applying an external load x/{(l — 
P{x))N). 

Let us first consider the breaking of fibers belonging to 
class 1 with threshold values uniformly lying between 
to cTi, this contribution £'i(A) is given by 



A(A) = 



A 



A-l 



A! 



dx 



1 



(l-a(x)) 



1 - ffi - .92 a{x) 
X exp({-a(a;) + In a(a;)}A) (11) 



where a{x) — x / {1 — gi — g2 — x) . The maximum contri- 
bution of this integral is at a{x) = 1 or a; = (1 — 51 — .g2)/2 
which exceeds ui, i.e., lies beyond the range of integra- 
tion and we can not employ the method of the saddle 
point integration. However, a{x) is a monotonically in- 
creasing function of x up to ai and we therefore get the 



Q 



O 




%A 



FIG. 2: Fig. (2a) corresponds to the ASD for easel of the 
table where a crossover from nonuniversal to the universal 
5/2 behavior is observed since 0-4 is close to 0.5. (2b) shows 
caseS of the table where no non-universal behavior is observed 
except for small values of A. ASD for a critical distribution 
(case4) is shown in Fig. 2(c) where there is a crossover to the 
"3/2" behavior for large A. The dotted line in Fig 2a and 2b 
has a slope -5/2 whereas that in Fig 2c has a slope -3/2. In 
2(d), we show Di{A), D2{A), D3(A) and the total ASD for 
the easel to compare their relative magnitudes. 



maximuni contribution when x reaches the upper limit of 
the integration cti. The leading behavior of Di{A) (with 
A! = exp(-A)A^V27rA), is given by 



Di(A) 



1 



2tt{1 - 51 - 52) 



A-'/^xteMa-^m)A),{12) 



where Xm = (Ti/{l-gi-g2-cri) = /i/(l-/i). Therefore, 
-Di(A) exhibits a non-universal decay with increasing A 
which is more rapid if /i ^0. 

Similarly, the leading contribution of the fibers belong- 
ing to the class 2 with threshold values ranging from 0-2 
to (73 and a{x) = x/{l — .92 ~ 2;) is found to be 



D2{A) 



1 



-5/2 A 



27r(l - 51 - 52) 



y;^exp((l-y™)A),(13) 



where ym = o'3/(l — <^4)- Therefore, D2{A) shows a sim- 
ilar non-universal behavior which survives even for rel- 
atively higher values of A if as approaches (T4 and also 
cr4 —^ 0.5. The significance of the above findings is ex- 
plained below. 

Let us now focus on the contribution from the fibers 
belonging to class 3. Using Eq. (10) with a{x) = x/{l—x) 



and /s == (1 - cr4)/(l - .91 - 92), we get 
1 /3 



D3iA) 



%/27rA 1-0-4 .L. X 



exp({- 



+ ln( )}A) (14) 



1 — a; "■ 1 — x 

Right hand side of Eq. (14) can be integrated to obtain. 



D3{A) 



h 



2^/2^A5/2(l-o•4) 



1 



(15) 



where Ac ~ (1/2 — ai)^'^, which diverges in the limit 
cr4 —> 0.5. Our observations are depicted in Fig. 2, where 
D(A), obtained by numerical integration and also by sim- 
ulation using the weakest fiber approach [1] , is plotted for 
different threshold distributions. It is to be noted that 
the contribution I?3(A) depends only on 0-4 and /s, an 
observation that leads to an interesting conclusion that 
the contribution to the total avalanche size distribution 
coming from the strongest class of fibers is identical for 
any number of discontinuities if the fraction of fibers as 
well as the range of the threshold strength of the final 
block is kept fixed (see Fig. 3). Equation (15) provides 
two limiting power-law behaviors given as 



DsiA) - A-3/2 for A < Ae 
- A-5/2 foj. A > Ac 



(16) 



Comparing Eqs. (12), (13) and (16), we observe that 
in the limit of small A, the non-power law contribu- 
tions from Di(A) and D2(A) dominates over the uni- 
versal "5/2" behavior only in the limit of tT4 -^ 0.5 when 
Z?3(A) ~ A~'^/2^ Otherwise, £'3(A) dominates over the 
non-universal contributions so that one observes a uni- 
versal behavior nearly for the entire range of A though 
the discontinuities in the distribution always exist. At a 
critical distribution however, there is a crossover to the 
universal behavior -D(A) ^ A^'^/^ foj- very large A fol- 
lowing a large region of non- universality (see Fig. 2c). 
This general behavior is valid for any number of discon- 
tinuities including n = 1 [l5|. 

Let us now concentrate on some interesting limiting sit- 
uations to investigate the role of two discontinuities: (i) 
If the fraction /i of fibers in class 1 is small, the contribu- 
tion Di(A) - [/i/(l - fi)]^ dies off rapidly However, if 
(73 is large and a^ approaches 0.5, there is a wide region of 
non-universality in (Z?(A) — A) behavior which is solely 
due to D2{A) which scales as [0-3/(1 — 0-4)]^. (ii) In the 
other limit, when /2 < /i (but /2 + /1 < 0.5, as required), 
Di{A) do dominate in the small A limit, but for large 
A once again it is the contribution of 1)2 (A), rather the 
larger value of 0-3 that leads to a prominent non-universal 
behavior, e.g., case 6 in table 1. We therefore conclude 
that the contribution from the weakest class of fibers is 
not significant in the large A limit and it is the higher 
value of Urn = 0-3/(1— 0-4) that causes the non-universality 
to survive up to relatively higher values of A (Fig. 4) 
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IV. CONCLUSIONS 



FIG. 3: Comparison of total avalanche size distribution for 
same value of fn+i and a^n when n — 1 and n = 2 where n is 
the number of discontinuities. Clearly, the two cases overlap 
in the large A region where the final block dominates. Here, 
/„+i=0.9 and cr2„=0.28 




lo^/ 



FIG. 4: The figure shows the comparison of D2{A) with 
D3{A) for case 1 (left panel) and case 6 (right panel). The 
figure shows that the crossover to universal behavior occurs 
at higher value of A if ym = 173/(1 — CT4) increases. 



In conclusion, we have studied a mixed fiber bundle 
with many discontinuities in the threshold distribution 
and GLS where threshold values of the fibers belonging 
to a particular class are uniformly distributed within the 
specified range. Our studies lead to the following conclu- 
sions for an arbitrary number of discontinuities: (i) The 
recursive dynamics studies point to the existence of a crit- 
ical distribution as defined in the text, (ii) There exists a 
non-universal, non- power law behavior in the avalanche 
size distribution for small A which becomes prominent 
only when a critical distribution is approached, otherwise 
it is masked by the universal behavior except for very 
small A. For asymptotically large A, however, there is 
always a crossover to the universal behavior with ^ = 5/2 
(or ^ = 3/2 at the critical distribution). The crossover 
occurs around A = Ac where the contribution from the 
strongest class of fibers D„+i(A) switches from ^ = 3/2 
to ^ = 5/2 behavior, (in) D„+i(A) is found to depend 
only on fn+i and a2n so the contribution of the fibers 
belonging the strongest class (i.e., the behavior of to- 
tal -D(A) in the limit of large A) remains identical for 
any number of discontinuities if fn+i and tT2n are kept 
fixed, (iv) We also show that if a2n-i increases and at 
the same time a2n ^ 0.5, the non-universality survives 
up to higher values of A. 
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